DESY 93-094 
HD-THEP-93-28 
July 1993 



m 
o\. 

Critical Exponents 
from the Effective Average Action 

< 

[ N. Tetradis 

Deutsches Elektronen- Synchrotron DESY, Gruppe Theorie, 
Notkestr. 85, 22603 Hamburg 



and 
C. Wetterich 



(N 
oo 

o 

Institut fur Theoretische Physik, Universitdt Heidelberg, 
Philosophenweg 16, 69120 Heidelberg 

CD 
> 

X 

.«?. 

Abstract 

We compute the critical behaviour of three-dimensional scalar theories using a 
new exact non-perturbative evolution equation. Our values for the critical exponents 
agree well with previous precision estimates. 



1. Introduction 



A new exact non-perturbative evolution equation for the scale dependence of an effective 
action has been proposed recently H . The effective average action Tk obtains by integrat- 
ing out only the quantum fluctuations with momenta q 2 > k 2 . Then the usual effective 
action (the generating functional of the 1PI Green functions) is recovered in the limit 
k — > where all quantum fluctuations are included. The effective average action can also 
be viewed as the action for averages of fields, similar in spirit to the block spin action 
|2|. |3| m lattice theories. The dependence of the effective average action on the scale k 
is described by the exact evolution equation. This equation has the form of a one-loop 
expression, differentiated with respect to k. The non-perturbative content, which makes 
this equation exact, consists of using the exact inverse propagator, as given by the second 
functional derivative of the effective average action T^, instead of the inverse propagator. 
Furthermore, an infrared cutoff suppressing the fluctuations with q 2 <C k 2 is added to the 
inverse propagator. In consequence, one obtains an equation where the scale dependence 
of Tfc is expressed in terms of its second functional derivative with respect to the fields 
(t = Ink) 

Ir^iTrf^ + i?,)- 1 !^}. (1.1) 

(Here is the effective infrared cutoff - details will be given in the next section.) The 
momentum integration implied by the trace is both ultraviolet and infrared finite in 
arbitrary dimensions. This makes our non-perturbative equation suitable for dealing with 
theories which are plagued by infrared problems in perturbation theory. (Scalar theories 
in less than four dimensions or at non-vanishing temperature near a second order phase 
transition and non-abelian gauge theories figure among the examples.) 

The most important new feature seems to us not the exactness of eq. (|1.1|) , but its 



simple and intuitive form which makes it suitable for approximate solutions. Exact equa- 
tions have been derived earlier, as for example the Schwinger-Dyson equations || or the 
"exact renormalization group equations" || ||. (There exists actually a formal relation 
between the latter and eq. ( |1.1|) 0-0. There should also be a formal relation to the 
Schwinger-Dyson equations - after all, exact relations describing the same Green functions 
should all be equivalent.) The evolution equation ( |1.1| ) is a functional differential equa- 
tion. It can be rewritten as an infinite system of ordinary coupled non-linear differential 
equations for the infinite number of couplings appearing in the most general form of Tk 
consistent with the symmetries 0. (For example, the general potential of a scalar field 
theory contains terms ~ <p 6 , ~ <f) 8 etc.) An exact solution of this infinite system seems 
almost impossible for an interacting theory. The r.h.s. of eq. ( |1 . 1|) involves the exact 
propagator which is as difficult to compute as Tk itself. The main point is that for many 
situations of interest the propagator is approximately known. Then one can employ the 
strategy to describe the propagator in terms of a finite number of parameters a.i which 
may depend on k. Typically, such parameters are masses, wave function renormalization 
constants etc. This approximation corresponds to an ansatz for the average action in 
terms of a finite number of undetermined functions ai(k). The evolution equation ( |1 . 1|) 
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is now truncated and turned into a finite number of differential equations 

(1-2) 

Such a finite system can be solved numerically or sometimes analytically. The sucess 
of our method ultimately depends on a good guess about the general form of the exact 
propagator, which in turn depends on the proper choice of degrees of freedom. (Sometimes 
it may be useful to include composite fields IjlQfl .) Fortunately, eq. ( |1.1| ) has a form 
very similar to the one-loop expression in standard perturbation theory. (It reproduces 
perturbation theory in its range of applicability. The loop expansion can be recovered 



from eq. ( |1 . 1| ) by a systematic expansion of the solution in terms of the small coupling.) 
This allows to build on the experience gained from the perturbative renormalization group 
analysis |4j], JTTJ] - [Tj| and provides useful checks for situations of overlap. Nevertheless, 



our equation is genuinely non-perturbative. Scalar theories at the critical temperature or 
in two and three dimensions are not accessible to perturbation theory and are described 



by our method with success [14], [L5] . 



For a scalar theory the inverse propagator is characterized by a mass matrix and a 
generalized kinetic term ~ Zq 2 . The mass matrix is directly related to derivatives of the 
scalar potential Uk which describes the non- derivative part of IV In general Z depends 
on fields and momenta and its detailed structure may be quite complicated. Nevertheless, 
this quantity plays essentially the role of a field and momentum dependent wave function 
renormalization. As a result its dependence on the various mass scales of the model is 
governed by the anomalous dimension. In particular, one expects a weak dependence on 
fields, momenta and the scale k if the anomalous dimension 77 is small. This is typically 
the case in three- and four-dimensional scalar theories. We will exploit this fact in the 
present paper. 

The aim of this paper is a precision test of our method in a context which is not 
accessible to perturbation theory. We consider the three-dimensional 0(iV)-symmetric 
scalar theory in the phase with spontaneous symmetry breaking. The presence of the 
Goldstone bosons leads to severe infrared divergences in perturbation theory and excludes 
a perturbative treatment in three dimensions. On the other hand, the model exhibits 
a second order phase transition. The critical exponents at the phase transitions have 



been computed by various other methods [12, [HJ with high precision. They serve as a 



valuable test. We shall see that even a truncation to only a few parameters gives good 
quantitative results for the critical exponents, typically with an accuracy at the percent 
level. In view of the fact that the momentum integration in eq. ( |1 . 1| ) is the same as for 
a one-loop calculation, and that no small coupling characterizes the interactions at the 
phase transition, this result seems quite remarkable. 

Since all our computations are done directly in three dimensions, we are not limited 
to the vicinity of the phase transition. In particular, we can employ our method in the 
case of the four- dimensional scalar theory at non-vanishing temperature and follow the 
transition from the effective four-dimensional behaviour for > T to an effective three- 
dimensional behaviour for k < T [[15] . The second order character of the high temperature 
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transition in 0(iV)-symmetric scalar theories has been established previously |T3|, together 
with a crude estimate of the critical exponents. The present work can be seen as a 
refinement of near the critical temperature where the three-dimensional behaviour 
dominates. A more accurate quantitative description of the temperature dependence of 
masses, expectation values and couplings away from the critical temperature, or the study 
of the full temperature dependent scalar potential will require to redo the computations 
of [|15|], including the additional couplings discussed in the present paper. 

This paper is organized as follows: We present the general formalism in section 2 
and derive in section 3 an exact evolution equation for the scalar potential formulated in 
terms of dimensionless couplings. No explicit mass scale appears in this equation, which is 
therefore well suited for a description of the scale independent physics at the second order 
phase transition. In section 4 we approximate this equation with a uniform (field and 
momentum independent) wave function renormalization. This approximation becomes 
exact in the large N limit. In section 5 we solve the exact evolution equation in the large 
N limit and establish consistency with earlier results [TO, llTj. In section 6 we discuss the 



general form of the effective scalar potential at the phase transition for arbitrary values 
of iV by taking the limit k — > for fixed values of the fields. In section 7 we exploit the 
knowledge that the anomalous dimension is small and sketch a truncation adapted to this 
property. The anomalous dimension is finally computed in section 8. Our results for the 
critical exponents are presented in section 9 and section 10 gives our conclusions. 



2. Exact evolution equation for the effective average 
action 

The scale dependent effective action T k , for a theory described by a classical action S, 
results from the effective integration of degrees of freedom with characteristic momenta 
larger than a given infrared cutoff k. The dependence of the effective action on the scale k 
is described by an exact evolution equation. In this section we summarize the formalism 
for the implementation of the above ideas in the case of a 0(iV)-symmetric theory of real 
scalar fields. For a detaided discussion, proof of the exactness of the evolution equation 
and modification of the conventions for the discussion of complex fields, we refer the 
reader to ref. 0. For our conventions see also the appendix A of ref. @. 

We consider a theory of N real scalar fields x a -, m d dimensions, with an O(N)- 
symmetric action S[x\- We specify the action together with some ultraviolet cutoff A, so 
that the theory is properly regulated. We add to the kinetic term a piece 

AS = ±Jd d qR k (q)x:(q)x a (q). (2.1) 

The function R k is employed in order to prevent the propagation of modes with charac- 
teristic momenta q 2 < k 2 . This can be achieved, for example, by the choice 

D , x Z k q 2 f^(q) 

Mq) = wife) ' (2 ' 2) 
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with 

/ fe 2 (g)=expl-^j. (2.3) 




As a result the inverse propagator for the action S + AS has a minimum ~ k 2 . The modes 
with q 2 ^> k 2 are unaffected, while the low frequency modes with q 2 <C k 2 are cut off: 

lim R k = Z k k 2 . (2.4) 

The quantity Z k in eq. ( |2.2| ) is an appropriate wave function renormalization whose 
precise definition will be given below. We emphasize at this point that the above form of 
R k is not unique and many alternative choices are possible. We subsequently introduce 
sources and define the generating functional for the connected Green functions for the 
action S + AS. Through a Legendre transformation we obtain the generating functional 
for the 1PI Green functions Tfc[0 a ], where (f) a is the expectation value of the field \ a m 
the presence of sources. The use of the modified propagator for the calculation of F k 
results in the effective integration of only the fluctuations with q 2 > k 2 . Finally, the scale 
dependent effective action is obtained by removing the infrared cutoff 



\jd d qRk{qW a {<l)nq)- (2.5) 
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For k equal to the ultraviolet cutoff A, T k becomes equal to the classical action S (no 
effective integration of modes takes place), while for k — > it tends towards the effective 
action T corresponding to S (all the modes are included). The interpolation of between 
the classical and the effective action makes it a very useful field theoretical tool. The 
means for practical calculations is provided by an exact evolution equation which describes 
the response of the scale dependent effective action to variations of the infrared cutoff 
(t = ln(*/A))i: 

|r fc [0] = ^Tr{(rt 2) [0] +R k )-^R k }. (2.6) 

(2) 

Here T k is the second functional derivative of the scale dependent effective action with 
respect to (p a . For real fields in momentum space Q, it reads 

< r ^'')=diko- (2 - 7) 

with 

n-q) = €(<?)• (2.8) 
Note also that, although not explicitly indicated for simplicity, R k in eq. ( |2.6|) stands for 

(R k ) a b (q,q') = R k (q)6 a b 6(q-q'). (2.9) 



* In order to avoid excessive formalism, we shall use the same notation for functions or 
operators in position space and their Fourier transforms, defined with the convention 4>{q) = 
(2ir)~2 J d d x exp(iq tl x >1 )4>(x) . 
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Tr 



°E/t&. (2-10) 



In momentum space, for a system whose volume f2 is taken to infinity at the end, the 
trace reads 

d d q 

and it is assumed to act on the diagonal part of the generalized matrix. The exact evo- 
lution equation gives the response of the scale dependent effective action to variations 

(2) 

of the scale k, through a one-loop expression involving the exact inverse propagator T k 
together with an infrared cutoff provided by R^- It has a simple graphical representation 
(fig. 1). Our non-perturbative exact evolution equation can be viewed as a partial dif- 
ferential equation for the infinitely many variables t and 4> a (q)- Its usefulness depends on 
the existence of appropriate truncations which permit its solution. 



Before presenting the formalism which leads to approximate solutions of eq. (|2.6|) , 
we briefly comment on the relation of the scale dependent effective action to the average 
action, which has been used for practical calculations in the past. The average action 
was introduced in ref. |7], [14] in order to describe the dynamics of averages of fields 
over volumes ~ k~ d . The implementation of an infrared cutoff k was naturally incorpo- 
rated through the averaging procedure. The dependence of the average action on k was 
computed through an one-loop renormalization group equation which corresponds to a 
truncation of eq. ( |2.6| ). It was shown in ref. || that a modification of the averaging 
procedure leads to an improved average action which is identical to the scale dependent 
effective action up to an explicitely known ultraviolet regulator in the improved average 
action. For this reason we call Tk the effective average action. Within the approximations 
used so far there is no difference for low momentum quantities between the originally 
proposed average action [0, 114] , the improved average action |§ and the effective av- 
erage action |l]. As a result, all previous calculations can be viewed as approximate 
solutions of eq. (|2.6| ) with some appropriate truncation. More specifically, these previous 
studies have described correctly the phase structure of the two- and three- dimensional 
scalar theories (including the Kosterlitz-Thouless phase transition) [L4|], the high tem- 
perature phase transition of the four-dimensional theory (with a reliable description of 
the effectively three-dimensional critical behaviour and a rough determination of critical 
exponents) fi~5| , |T7l , and the approach to convexity for the effective potential in the phase 
with spontaneous symmetry breaking |L8 



As we have already mentioned, for the solution of eq. ( |2.6|) one has to develop an 
efficient truncation scheme. We consider an effective action of the form 

T k = J d d x{u k (p) + ^dy a Z k (p } -d v d v )d^ a + ^pY k (p, -d v d v )d»p), (2.11) 

where p = \<p a ( t )a an d appropriate normal ordering is assumed for the derivative operators 
]![]. In order to turn the evolution equation for the effective average action into equa- 
tions for Uki Zk and Y^, we have to evaluate the trace in eq. ( |2.6|) for properly chosen 
background field configurations. For the evolution equation for Uj. we have to expand 
around a constant field configuration. This calculation is carried out in section 3 and the 
study of the resulting equation occupies sections 4-7. The evolution of the wave function 
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renormalizations, which leads to the determination of the anomalous dimensions, requires 
an expansion around a background with a small momentum dependence. This calculation 
is presented in section 8 and appendix B. The evolution equations for U k , Z k and Y k are 
partial differential equations for independent variables t and p. In most cases their study 
is possible when they are turned into an (infinite) set of coupled ordinary differential 
equations for independent variable t. This is achieved by Taylor expanding U k , Z k and Y k 
around some value of p. Since in this work we are interested in the vacuum structure of 
the theory, we shall use an expansion around the k dependent minimum po(k) of U k . In 
the limit k — > this minimum specifies the vacuum of the theory, and U = Uq, Z = Zq and 
Y = Y and their p-derivatives at the minimum give the renormalized masses, couplings 
and wave function renormalizations. 

In a first approximation we can neglect the q 2 dependence of Z k (p,q 2 ). We can now 
define the quantity Z k appearing in eq. ( |2.2| ) as 

Z k = Z k {p (k)). (2.12) 

For studies which concentrate on the minimum of U k , the above definition permits the 
combination of the leading kinetic contri 
inverse propagator (for massless fields) 



combination of the leading kinetic contribution to T k and R k in eq. (|2.6|), into an effective 



P(q 2 ) = Z k q 2 + R k = 1 ^ ■ (2.13) 

For q 2 ^> k 2 the inverse "average" propagator P{q) approaches the standard inverse 
propagator Z k q 2 exponentially fast, whereas for q 2 <C k 2 the infrared cutoff guarantees 
P(0) = Z k k 2 . Throughout most of this paper, we shall work with the approximation of 
q 2 independent Z k (p) and use the definition of eq. (|2.12| ). When the q 2 dependence of 



Z k (p, q 2 ) is considered, eq. ( |2.12|) has to be replaced by an appropriate generalization, for 
example 

Z k = Z k (p {k),q 2 = 0). (2.14) 

Having summarized the basic formalism, we obtain in the next section the exact 
evolution equation for the effective average potential U k . 



3. Exact evolution equation for the effective average 
potential and its scaling form 



The evolution equation for U k (p) can be obtained by calculating the trace in eq. 
for small field fluctuations around a constant background configuration <fi a {q = 0) = 
<Mai> P — \<t> 2 - We define 

(j) a (q) = (j)5 al + Xa(q) (3.1 

and find 

(rf <?') = [(Z k (p, q 2 )5 ab + pY k (p, q 2 )5 al 5 bl ) q 2 + M 2 ab ] 5(q - q'\ 



3.2) 
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where the mass matrix M 2 b reads 



M 2 ab = U' k (p)6 ab + 2pU'M6 al 6 bl . (3.3) 
Here primes denote derivatives with respect to p. Eq. (|2.6j) now reads 



d TT . . 1 f d d q d „ , . (N - 1 1 . 

m UM = 2J WT d dt Rkiq) \-m^ + w^ (3 - [) 



with R k (q) given by eqs. O - Q, 

M (p, g 2 ) = Z fc (p, g 2 )g 2 + i2 fc (g) + C^(p) 

M 1 (p,g 2 )= Z fc (p,g 2 )g 2 + i? fc (g) + ^(p) + 2p^(p) (3.5) 

and 

Z k (p, q 2 ) = Z k (p, q 2 ) + pY k (p, q 2 ). (3.6) 

For p different from zero it is easy to recognize the first term in eq. ( |3.4| ) as the contribution 
from the N — 1 Goldstone bosons (U' vanishes at the minimum). The second term is 
then related to the radial mode. The effective average potential U k has a k dependent 
minimum po(k). We introduce the wave function renormalizations for the Goldstone and 
radial modes (in agreement with eq. ([2.14 )) 

Z k = Z k (p (k),q 2 = 0) 
Y k = Y k (p (k),q 2 = 0) 

Z k = Z k (p (k), q 2 = 0) = Z k + p (k)Y k (3.7) 

When the q 2 dependence of Z k (p,q 2 ), Z k (p,q 2 ), Y k (p,q 2 ) is neglected the wave function 
renormalizations are defined according to eq. ( |2.12| ) and read 

Z k = Z k (p (k)) 
Y k = Y k (p (k)) 

Z k = Z k (po(k)) = Z k + p (k)Y k . (3.8) 
The k dependence of these functions is given by the anomalous dimensions 

m= -j t ^z k ). (3.9) 

At this stage the r.h.s. of eq. (|3.4|) depends explicitely on the scale k once the momentum 
integration is performed. The aim of this paper is a study of the critical behaviour at 
and near the second order phase transition. At the phase transition one expects a scaling 
behaviour of the effective average action Y k and, in particular, of the effective average 
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potential U k . This holds provided approximate renormalized dimensionless couplings are 
used. By a proper choice of variables the evolution equation ( |3.4j) should therefore be 
cast into a form where the scale k no longer appears explicitely. In this formulation it 
will become easy to study the fixed point behaviour. We perform the transformation of 
eq. (|3.4|) into its scaling form in several steps. First we introduce the variable x = q 2 and 
define the dimensionless functions 

Rk{x) 



r k {x) 



s k (x) 



Z k x 

]_d_ ( R k (x) 
Z k dt 



x 



d d 

— r k (x) + —(\nZ k )r k (x) 



d 

-2x—r k (x) - T)r k (x)(3.10) 
ox 



Here R k (x)/Z k x is a dimensionless function of the ratio xj k 2 (a property of R k independent 
of the specific choice (|2.2|)). Eq. (|3.4|) can then be written in the form 



d_ 

dt 



U k (p) =v d dxx*s k 



N 



+ 



with 



2 d+1 vrfrl- 



M /Z k MJZ, 



d 



The ratios Mo/Z k and M\jZ k read 

M (p, x)/Z k = z k (p, x)x + r k (x)x + Z k l U' k (p) 

M 1 (p,x)/Z k = zz k (p,x)x + r k (x)x + Z k l U' k (p) + 2Z k 1 pU k (p), 

where 

z k (p,x) = Z k (p,x)/Z k 
z k (p,x) = Z k (p,x)/Z k 

and 

z(k) = Z k /Z k . 

This suggests the use of a dimensionless renormalized variable 

P 

In terms of p one obtains 



Z k k 2 ' d p. 



dU k (p) 



dt 



(d-2 + V )p^M + v d k~ 2 



dp 



dxx 2 s k 



where we have used the identity 



dU k 



dt 



dU k 



Dt 

dUk 
dt 



+ 



dU, 



dp 



dp 
dt 



dU k 



dt 



+ P 



dU k 



dp 



N - 1 
m 



dlnp 



dlnp 



1 

mi 



dlnp 



dt 



+ (d-2 + r])p 



dUk 
dp 



(3.11) 
(3.12) 



(3.13) 



(3.14) 
(3.15) 

(3.16) 
(3.17) 



(3.18) 
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The dimensionless ratios mo and m\ are given by 

m (p,x) = = ^ [z k (p,x) +r k (x)} + k~ d ° l r[r ' 



Z k k 2 k 2 1 vr ' ' K n dp 



™i{p, x) = = [z2k(p, x) + r k (x)] + k d 



dUk(p) | 2 ~ d 2 U k (p) 
dp dp 



n2 



[3.19) 



We finally switch to a dimensionless potential 

u k (p) = k- d U k (p) (3.20) 

and use a dimensionless momentum variable 

y = x/k 2 . (3.21) 

We adopt the convention that the prime on u k (p) denotes a derivative with respect to 
p at fixed t, while the prime on U k (p) denotes a derivative with respect to p at fixed t. 
At this point it is useful to fix the notation which we shall use in the next sections. We 
denote the minimum of U k (p) with po{k) and the minimum of u k (p) with K{k) = po(k). 
For po{k) 7^ (spontaneously broken regime) we have U k (po) = 0, u' k (K) = 0. We define 
UZ(po) = X(k), U { k n \p Q ) = U n (k) and = X(k), u { k n \ K ) = u n (k). For p (k) = 

(symmetric regime) we define U k (0) = rh 2 (k), u' k (0) = m 2 (k)/k 2 = m 2 (k) and as above 
for the higher derivatives. We point out the relations 

k = Z k k 2 ~ d p 

~ 2 ry-l r -2 - 2 r-2 2 

m = Z k k m = k m 
X= Z k 2 k d ~ A \ 

Z k n k {n - 1)d - 2n U n . (3.22) 



u, 



We can now write the evolution equation in a scale independent form 

d t°° a ( N — 1 1 \ 

Ol u k(p) = -du k {p) + (d-2 + rj)pu' k {p) + v d ^ dyy*s(y) [~^- + —) , (3-23) 

with 

m (p,y)= y[z k (p,y) +r(y)}+u' k (p) 

mi (p, y)= y [zz k (p, y) + r(y)) + u' k {p) + 1~pu!' k {p). (3.24) 

Indeed, r k = r(y) only depends on y, whereas s k = s(y,rj) = —2y^r[y) — i]r(y) depends 
on y and rj. We have, therefore, dropped the label k for these functions. 

Eq. ( |3.23|) is the scaling form of eq. ( |3.4|) we were looking for. The behaviour exactly 



at the second order phase transition should be given by a /c-independent solution of this 
equation. The scaling solution of (|3.23| ) obtains therefore for du k /dt = 0. Since no explicit 
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k dependence is present on the r.h.s. of ( |3.23D we conclude that all the k dependent 
quantities in ( |3.23|) obtain constant (fixed point) values. In particular r](k), fj(k) and 
z(k) take values 77*, 77* and for this solution. The definitions ( |3.9|) , ( |3.15|) then imply 
77* = 77*. Similarly w*(p) and z*(p,y), z*(p, y) become functions which are independent 
of k. The nature of the fixed point becomes apparent through the observation that the 
minimum of the potential appears for a constant value 

K(k) = «*. (3.25) 

As a result, for the fixed point solution, the minimum of the effective potential is given 
by 

po = limpo(fc) ~ lim k d ~ 2+ri * = 0. (3.26) 

We conclude that the fixed point of the evolution equation indeed corresponds to the phase 
transition between the spontaneously broken and the symmetric regime. The scaling form 
of the potential w*(p) is determined by the equation 

du+{p) - (d - 2 + V *)P<(P) + 2v d r dyyi (y^- + ^r(y)) ^ + — ) = 0, 

Jo \ ay 2 ) \ m * m u J 

(3.27) 

with 

mo*(p,y)= y[z*(p,y) + r (y)} + u '*(p) 

mu(p, y)= y [z*z*{p, y) + r{y)\ + u'JJp) + 2pu"(p). (3.28) 

If the functions z*(p, y), z*(p, y) and the constants 77*, 5* are known and the ^/-integration 
is performed, eqs. ( |3.27|) , (|3.28|) are reduced to a non-linear second order differential 
equation for u*(p). 

We point out that eqs. ( |3.23|) , ( |3.27| ) are exact, since eq. (|2.11|) contains the most 



(2) 

general terms which contribute to T k when evaluated at a constant background field. 
We have, therefore, obtained an exact non-perturbative evolution equation which does 
not involve any explicit mass scale. 



4. Approximation with uniform wave function renor- 
malization 



The solution w*(p) of the scaling equation ( |3.27p together with the evolution equation 



( |3.23p contains all the important information about the behaviour at and near the phase 
transition. Indeed, eq. ( |3.23| ) gives the functional dependence of the various /5-functions 
on the couplings at the phase transition. (In fact, u*(p) describes infinitely many couplings 
and eq. (|3.23|) specifies infinitely many /5-functions.) This allows to compute derivatives 
of the /5-functions with respect to the couplings. From there one may compute the matrix 
of anomalous dimensions and infer the critical exponents at the phase transition. (For 
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example, the critical exponent v is related to the anomalous mass dimension.) Even 
though eqs. (|3.23| ) and (|3.27p are exact the r.h.s. of eq. ( |3.27|) still involves the unknown 



constants 77*, z* and functions z*(p,y), z*(p, y). If these are given one may attempt to 
compute the y-integral in eq. ( |3.27D numerically. The result is a function depending on 
p, u'^(p) and u"(p). The remaining non-linear second order differential equation can also 
be treated with numerical methods. This program requires some effort, but the result - 
exact values of the critical exponents - would be rewarding. 

Unfortunately, exact information on 77* , 2*, z*(/5, y) and £*(/5, y) is not available and the 
problem seems only postponed. However, the situation is not so hopeless, at least for the 
three-dimensional models for which 77* is typically a small quantity and z^(p,y), z*(p,y) 
are expected to depend only weakly on p and y. Then, a relatively crude approximation for 
T[* , Z* , Z* (p, y), z*(p, y), which may be obtained by truncating the exact evolution equation 
for these quantities or by alternative methods, should already give excellent results. We 
shall demonstrate the power of our method by making the very simple approximation of 
uniform wave function renormalization, i.e. 

z k (p,y) = z k (p,y) = 1- (4.1) 

This corresponds to neglecting the field and momentum dependence of the wave function 
renormalization. We observe that the dominant scale dependence of the wave function 
renormalization - the anomalous dimension - is always included by the fc-dependence of 
Z k and Z k (the definition (|3.8|) applies). For non-zero values of the infrared cutoff k 



the more general wave function renormalization function Z k (p,q 2 ) should be analytic at 
q 2 = even at the phase transition. This is a consequence of the presence of the infrared 
cutoff. One expects a very mild momentum dependence of Z k (p , q 2 ) for q 2 between zero 
and k 2 whereas for q 2 ^> k 2 the anomalous dimension should determine the behaviour 
near the phase transition (Z k (p ,q 2 ) ~ (g 2 ) _ a). We note, however, that all momentum 
integrals are dominated by q 2 ~ k 2 . Similarly, an expansion of Z k (p,q 2 ) in the field p 
should be analytic at po(k) as long as k is strictly positive. (See section 6 for more details.) 
Our approximation ( |4.1| ) should, therefore, correctly describe all qualitative features at 
the phase transition and be a good quantitative approximation. Of course, if a high 
precision calculation is attempted one will have to go beyond the truncation ( |4.1| ). In the 
approximation ( fl.lp it is useful to define the integrals 

L d (w, 7),z)= 2 dyy^ +1 ^^- [y(z + r(y)) + w]' 1 

+ r) dyy*r(y)[y(z + r(y)) + w} 



L d n (w, 77, z) = 2n [°° dyy^ 1 ^- [y(z + r(y)) + w]~ {n+1) 
Jo ay 

/"CO 1 , . 

+ nrj dyy^r(y) [y(z + r(y)) + w]~^ n+ ' for n > 1, (4.2) 

Jo 
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which obey 

d ~ ~ 
—L d (w,rj,z) = -Lf(w,r},z) 

d ~ ~ 

—L d (w,rj,z) = -nL d n+1 (w,7],z) forn>l. (4.3) 
Similarly we introduce 

L d n (w, V ) = L d n (w, V ,l). (4.4) 

In the following, for notational simplicity, we shall not write down explicitly the depen- 
dence of L d and on 77 and z. Also, for the practical calculations in section 9, we 
shall omit the contribution to L^, L d proportional to r] (we shall study the critical three- 
dimensional theory, for which 77* ~ — 0.04) and we shall set z = 1 (which is again a 
good approximation near the fixed point). With these simplifications we need only the 
detailed form of L d (w) = L d (w,0). For the choice Q2.2p - ( |2.3j ) for Rk(q) these integrals 
are discussed in appendix A. We emphasize, however, that all our formulae remain valid 
for 77 7^ and z ^ 1, even though the dependence on these parameters is suppressed in 
the notation. 

In terms of L d , L d eq. ( |3.23|) reads 

_ = - du + id - 2 + V )pu' - v d (N - l)L d {v!) - v d L d {u' + 2pu"), (4.5) 

where, for simplicity, we drop the subscript k of Uk from now on. The evolution equations 
for the derivatives of u are easily obtained by differentiating (4.5) with respect to p: 

du' 



dt 
du" 

In 



(-2 + 77)14' + {d - 2 + v )pu" 

+ v d {N - l)u"L d {u') + v d {3u" + 2pu"')L d 1 (u' + 2~pu") (4.6) 
(d - 4 + 2i])u + (d - 2 + r])pu" 
-v d {N-l){[u"fL d {u')-u'"L d {u')} 

- v d {(3u" + 2pu'") 2 L d (u' + 2pu") - (5m w + 2pu {i) )L d l (u' + 2pu")). (4.7) 

In the spontaneously broken regime, at the minimum n(k) we have u'(k) = 0, A = u"(k) 
and u n = u^ n '(K). The running of the minimum is obtained by taking a total i-derivative 
of the condition u'(k) = 

^ = A = -KM]" 1 ^ 



p=K 



-(d-2 + v )K- v d (N - l)L d {0) -v d (s + L?(2Ak). (4.8) 

Eq. (|4.5|) remains a complicated partial differential equation for two variables t and p, 
even if L d , Lq take a simple functional form. If we consider u in the vicinity of some fixed 
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value of p, e.g. p = k, we can parametrize u by its derivatives at this fixed value, e.g. A, 
u n . The evolution equation is then expressed as an infinite system of coupled differential 
equations for the couplings X,u n . In particular, the /3-function for A reads 



dX du" 

-17 — PX 



dn 

p=K 



dt dt 

(d - 4 + 2 V )X - v d {N - 1)A 2 L^(0) 



/ 2km 2 \ ~ 
- v d (3X + 2ku 3 ) 2 L%(2Xk) + v d I 2u 3 + 2ku 4 ^ J L?(2Ak) (4.9) 

and it involves the higher derivatives u 3 and u 4 . In turn, the /3-function for w 3 contains 
m 4 and u 5 and so on (see section 7 for the details). It should be also pointed out that 
the functions L^, automatically introduce a "threshold" behaviour in the /3-functions, 
which leads to the vanishing of the contributions from the massive radial mode at scales 
k much smaller than its running mass. 

As a starting point for the understanding of the solutions of ( |4.5| ) we need first the 
behaviour of the scaling solution for du/dt = 0. The scaling equation 

- du* + (d - 2 + 7j,)K - v d (N - l)L d (u'J - v d L d (u'* + 2p<) = (4.10) 
is a non-linear differential equation for p, which remains difficult to solve analytically. 



Numerical work for its solution is in progress [T^j. For the present paper we rather 
concentrate on analytical aspects and present, in the next two sections, results for limiting 
cases. 



5. Solution of the evolution equation for the effective 
average potential for N —> oo 

The main results of this work will be obtained through use of the evolution equation 
( |4.5|) for the study of the three-dimensional theory. This will require a certain amount of 
numerical work. It is useful, therefore, to have analytical results in some limiting case, in 
order to obtain conceptual understanding and verification of the numerical solution. This 
is the purpose of this section in which we concentrate on the three-dimensional theory in 
the large N limit [pOfl . For iV — ► oo the evolution equation ( |4.5|) simplifies considerably. 
The first simplification obtains from the dominance of the Goldstone contributions in the 
evolution equation. The second one results from the vanishing of the anomalous dimension 



rj. (This is a known result [|T^] which will be reproduced in section 9). This means that we 
can set Zj, = 1 and omit all terms involving rj in our expressions. It is more convenient to 
study the equation (|4.6|) for vl . The potential can be obtained at the end through simple 



integration. In the large N limit eq. ( |4.6|) takes the form (d = 3) 



^- = -2v! + pu" + VzNu"L\{u'). (5.1) 
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For the choice of Rk of eqs. ( |2.2| ), ( |2.3[ ) the function L\{u') does not have a simple 
analytical form. We use the notation of eq. (|A.5|) and and make, in this section, the 
approximation (with l\ a constant of order unity) 

L\(u')= -2l\s\(u') 

s\(u') = {l+u'Y 2 . (5.2) 

Notice that the above choice preserves the correct behaviour of L\(u') as discussed in 
appendix A: monotonicity, existence of a pole at u' = —1, and asymptotic dependence 
on u' for u! — > oo. We emphasize that the choice of Rk is arbitrary within some general 
conditions It may be possible to find an Rk such that eqs. ( |5.2|) become exact, with 
a suitably chosen value of If. (See also footnote on page 17.) 
Eq. ( |5.1|) can now be written as 

dv! „dv! NC du' n , n 

ln- p W + aTWW + 2u °> (5 - 3) 

where C = 2vslf = If/Air 2 . We have also kept explicitly the factor of JV for consistency 
checks through the N dependence of the couplings. The most general solution of the 
partial differential equation ( |5.3|) is given by the relations 



p NC NC VW 3 Ar ^ / 1 \ „ f / »\ , / , s 

' 1 -AC arctan —= = F (u'e 2t ) for u' > (5.4) 



w VW 2 1 + u' 2 \Vu' , 

—£= ^=H AC In -== }=F(u'e 2t ) for u' < 0.(5.5) 

V^u 1 2 1 + u' 4 \i + ^^J \ J v > 

The function F is undetermined until boundary conditions are specified. For t = 
(k — A), £7fc coincides with the classical potential which we choose 

^A(p) = ^) = -^ + yP 2 (5-6) 
The boundary condition, therefore, reads 

u'(p,t = 0) = -p J 2 A + \ A p, (5.7) 

where 

A = |. A A = ^. (5.8) 
This uniquely specifies F and we obtain 



p NC NC VW 3 Ar ^ 

H — NC arctan 



u' VW 2 1 + u' 2 \Ju' 



We 1 (p\ \ 1 NC y/v/e* 3 ^ / 1 \ 

+ - JVC — — - , nj + -JVC arctan for u' > 



Aa VA a yfufe 1 2 l + w'e 2i 2 Vv^e* 
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NC NC v / -m 7 3 Ar ^, /l-\/^F 

-NC In 



^+^-^^ + ^^4-?^^-^ for „'«,. 

(5.10) 



A A \ A a ' / v^e* 2 l+w'e 2 * 4 \i + v /_ M - e 



Even though the approximation of eq. (|5.2j ) does not permit quantitative accuracy in 
all respects, eqs. (|5.9|) , ( |5.10|) contain all the qualitative information for the non-trivial 
behaviour of the three-dimensional theory. There is a critical value for the minimum of 
the classical potential 

«(A) = ^ = n cr = NC, (5.11) 
Aa 

for which a scale invariant (fixed point) solution is approached in the limit t — ► — oo. This 
solution is given by the relations 



p NC NC \FuL t 3 1 3tt ^ . , . _ 

1 x + -NC arctan -— = — NC for <> 0(5.12) 



K vX 2 l + < 2 Vv^Z 4 



^ ^ + - ?*Cln f 1 -^ = for < < 0.(5.13) 



2 l + < 4 + 

On the other hand, the general solution of the scaling equation, which results from eq. 
by setting du'/dt = 0, is given by 



NC N ° ^ ' 3 NC arctan = c + for < > (5.14) 



m' v^i 2 l + < 2 



^ + - ?^ln fi^l = c_ for < < 0, (5.15) 



V=< 2 l + < 4 ^i + v ^ 

with c+, c_ arbitrary constants. It is remarkable that there is only one solution of the 
scaling equation which is continuous and finite for all finite p. This is exactly given by 
eqs. ( 5.12 ), (|5.13j ). This criterion specifies more generally the selection of the two free 
integration constants which are necessarily present for the general solution of eq. ( |4.10 ) . 
Eqs. (|5.12|) , ( |5.13|) describe a potential u which has a minimum at a constant value 

K {k) = k* = NC, (5.16) 

and, according to ( |3.26|) , corresponds to the phase transition between the spontaneously 
broken and the symmetric regime. (The values for n cr and coincide, but this is acci- 
dental.) For the second p-derivative of u at the minimum A = u"{k) we find 

A(A0 = ^ = (5.17) 
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and similar fixed point values for the higher derivarives of u. For 1 <C p/NC <C (37r/4)e * 
the rescaled potential u has the form 



Notice that the region of validity of eq. (|5.18|) extends to infinite p for t — > — oo. From eq. 



( [5.18 ) with t — > —oo(k — > 0) we obtain for the effective potential at the phase transition 

and a critical exponent 5 = 5, in agreement with standard large N results [M. (For the 
definition of 5 and detailed discussion for non-zero rj see section 6.) 

Through eqs. ( |5.9| ), ( |5.10| ) we can also study solutions which deviate slightly from the 
scale invariant one. For this purpose we define a classical potential with a minimum 

k(A) = n cr + N6k a , (5.20) 

with \8k\\ <C 1. We find for the minimum of the potential 

K {k) = + N5K h e~\ (5.21) 



and for A 



M 

Eq. (|5.21|) indicates that the minimum of the potential stays close to the fixed point value 
given by ( |5.16|) , for a very long "time" \t\ < — In \ Sk\\. For \t\ > — In \Sk\\ it deviates 



X(k) = K (5.22) 

1 + ( ^ - 1) e* 



from the fixed point, either towards the phase with spontaneous symmetry breaking (for 
5ka > 0), or the symmetric one (for 5k,\ < 0). Eq. ( |5.22| ) implies an attractive fixed 
point for A, with a value given by eq. ( |5.17| ). Similarly the higher derivatives are attracted 
to their fixed point values. The full phase diagram corresponds to a second order phase 
transition. For 5k\ > the system ends up in the phase with spontaneous symmetry 
breaking, with an order parameter given by 

p = limpo(^) — limkn(k) = NSk^A. (5.23) 
This leads to a critical exponent 



(3 = lim L v j = 5 ^ ^ 5>24 ^ 



Sk a ^o d [In 5k^] 

again in agreement with standard large N results [|T2 1 . (For a detailed discussion of 



the exponents j3, is, ( see section 9.) In this phase the renormalized quartic coupling 
approaches zero linearly with k 

\ R = lim k\{k) = lim kK = 0. (5.25) 
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The fluctuations of the Goldstone bosons lead to an infrared free theory in the phase with 
spontaneous symmetry breaking. For 5k a < 0, n(k) becomes zero at a scale 



In 



N\6k a \ 



[5.26) 



and the system ends up in the symmetric regime (po — 0). From (|5.9| ), in the limit 

t -> -oo, with u r ,u",u'" -> oo, so that u'e 2t ~ \Sk a \ 2 , u"e l ~ \6k a \/N, u'" ~ 1/N 2 , we 
find 

4 



U(p) = U (p) 



3nC 



2 r i I 

[5k a ] 2 A 2 P - -5n A Ap 2 + — 



P 



[5.27) 



in qualitative agreement with previous studies [12], |T^| . In particular, we emphasize that 
the critical behaviour of the three-dimensional theory coincides with the behaviour of the 
four- dimensional theory at non- vanishing temperature near the critical temperature T cr . 
Eq. (|5.27| ) yields for both the symmetric phase and the phase with spontaneous symmetry 
breaking (p > p min ) 



U'( P 4 




provided we identify 



5k a A 



J rp2 rp2 



[5.28) 



[5.29) 



24 T 

and use the four-dimensional normalization for the field p^ = pT. Up to a possible 
difference [j] in the overall normalization of U' this is exactly the high temperature result 



of ref. fL7 |. One obtains the following well known values for the critical exponents z/, ( in 
the symmetric phase 



d [In m R ) 

lim : : 

Sk a ^o d [In \Sk a \] 

Um d[\n\ R ] 
Sk a ^o d [In \Sk a \] 



d 



lim 

<5«a^o d[ln|5«; A |] 

d[lnU"(0)] 
hm — F — — — — 

Sk a -*o d [In \6ka\] 



1. 



t inn — t — — — — ma — : , — --- = I (5.30) 

(5.31) 

«b < 5 - 32 > 

a result connected to the resolution of the problem of infrared divergences of perturbation 
theory 



Finally we find 



lim — — 



Another approach for the study of eq. (5^) would make use of a parametrization of 



u in terms of its minimum k and the derivatives at the minimum A, 



u, 



as discussed at 



t Eq. (3.15) in ref. jlTj obtains for i f = 2/3. It differs from eq. ( |5.28 ) in the present paper by an 
overall fact or if we insert the value ( A. 4 ), since If = This difference is du e to th e ap proximation 

of eq . (5.2). Note that it may be possible to alter the definition of Rk in eqs. ( |2.2| ), (2.3) so that eq. 
(5.2) becomes exact. Then If should come out to be 2/3 with this special choice. 
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the end of the previous section. For the approximation of eq. 
limit, eqs. flOft , Q read 



"(ft 

dA 



- k + NC 



- A + 2NCX . 



j^), and in the large iV 

(5.33) 
(5.34) 



The solution of the above equations, for a classical potential of the form of eq. ( |5.6| ), 
is given by eqs. ( |5.21| ), (|5.22| ). Notice, however, that the disappearance of the higher 
derivatives of the potential in eqs. (|5.33| ), (|5.34j) is a very convenient simplification 
appearing only for N — > oo. In the next sections, where we shall study the full eqs. ( |4.8| ), 
( |4.9| ) for small values of N, more elaborate approximation schemes will be necessary. 

As a final comment we point out that, for a theory with spontaneous symmetry break- 
ing , we can use eq. ( |5.10| ) in order to study the "inner" part of the potential. In particular, 
for p = and t — > — oo eq. ( |5.10| ) predicts a potential u which asymptotically behaves as 



lim 



«'(0) 



1. 



This leads to an effective average potential Uk which becomes convex with 

lim £7^(0) = -A: 2 , 



^5.35) 



^5.36) 



in agreement with the detailed study of ref. . 



6. Scaling solution for k — > at fixed p 



We now turn back to the finite values of N. The scaling equation ( |4.1Q|) takes a partic- 



ularly simple form for p — > oo which we will study in this section. We observe that the 
integrals Lq(w), Lq(w) rapidly decrease for w ^> 1 (see appendix A). For large values of 
p we, therefore, obtain the asymptotic equation 

-du Sf + (d-2 + ri*)pu'^ = 0. (6.1) 

This has the simple general solution 

u*(p) = cp T (6.2) 

with 

r=— — (6.3) 

and c an integration constant. 

The limit k — > for fixed p corresponds to p — > oo for d — 2 + 77* > (see eq. (|3.16Q ). 
In this limit Uk becomes the effective potential Uq = U. As a result eq. ( |6.2|) gives the 
scaling form of U(p), namely 

U^p) = cZlp\ (6.4) 
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Here the constant Z c appears in the scaling form of the wave function renormalization 



(6.5) 



The form of the potential (|6.4|) is directly related to the critical exponent 8 characterizing 



the relation between the magnetization <fi and the magnetic field B = Wr = <fi¥¥- 



■8U 



dp 



B 



~ (7) ~ 



(6.6) 



One obtains the well known scaling law 



5 



d + 2 — ri* 
d — 2 + rj* 



(6.7) 



At this point the reader may be somewhat surprised that the simple approximation 
automatically produces the correct relation between 5 and 77*. One would expect 
that for k d ~ 2 much smaller than Z k p the wave function renormalization Z k (p) becomes 
independent of k, since no low mass modes are present. The identification of Z k (p) with 
Z k as defined in eq. ( |3.8|) becomes then problematic. Indeed, the presence of two different 
infrared cutoffs: k 2 and 



m 



Z k \p) 



dU 
dp 



(6.8) 



suggests a qualitative behaviour 



Z k {p) ~ {k 2 + Cm 2 



2 



(6.9) 



This implies that Z k (p) and Z k = Z k (p (k)) differ substantially for p ^> 1. Nevertheless, 
the asymptotic form of the scaling equation, given by eq. ( |6.1| ), remains valid. This can 
be seen by introducing a different scaling variable 



p = Z k {p)k 2 - d p 



(6.10) 



which absorbs the dominant effects of the wave function renormalization for large p as 
well. Using the scaling behaviour given by eqs. (|6.4|), (|6.9|) for large p we find 



_ rf* (t — 1) r}* 

Z k {p) ~ m~ Th ~ p 2 -'<* ~ p d ~ 2 +r>* 



p ~ k a pd-2+ n * 



(6.11) 
(6.12) 



We replace eq. ( |3.18j ) by 



dU k 



Of 



dU k 



Of 

du k 



+ 



dU k 



dp 
dt 



dt 



dp 

+ (d-2)p 



dU k 



p 

dU k 



dt 



+ p 



dU k 



dp 



dlnp 



dhip 



dlnp 



dt 



dp 



(6.13) 
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The asymptotic solution obtains in a similar way as before 



(6.14) 



and coincides with eq. ( |6.4| ). We conclude that, in a leading approximation, the scaling 
equation ( |6.1|) remains valid for large p despite the difference between Zk(p) and Zk- 
We finally observe that corrections to the asymptotic solution can be obtained from 



eq. ( 4.10 ) by factoring out the asymptotic behaviour 

u{p) = u{p)p~ T . 



[6.15) 



This leads to 



(d - 2 + V *)pK-v d (N - l)p- T L d (u'J T + 



" ~T— 1 

TU^p 



-v d p- T L d (2KP T+1 + (1 + 4t)u'J t + r(2r - I 



u*p T 1 



0. (6.16) 



The integrals Lq, Lq can be expanded for large values of the argument w. (For the 
asymptotic behaviour of Lq for rj = see eq. (|A.6|) in appendix A.) In this way one can 
compute the leading corrections to eq. (|6.2| ) for large p. We postpone this discussion for 
a future publication J19 . 



7. The T] expansion 

In three dimensions, where 

r=— ^— (7.1) 

1 + 77* 

we observe that, according to eq. (|6.2| ), the third derivative of «*(/5) goes to zero for 
p — ► oo since i]* is positive 

u'l'(p) = cr{r - l)(r - 2)p~ — . (7.2) 



This also holds for the higher derivatives of u*. On the other hand we learn from eq. 
that U"'(p — > 0) diverges at the phase transition. This shows that the </> 6 coupling plays an 
important role at and near the phase transition. Since the fourth and higher derivatives 
of U* also diverge as p — >• one may suspect that all these higher derivatives have to be 
included for a quantitative estimate of the properties near the phase transition. This would 



require a solution of eq. ( |4.5| ) without the possibility of truncating the corresponding 
system of infinitely many differential equations for (see the discussion at the end of 
section 4) by setting = for n higher than a certain integer. 

Fortunately, the situation is much better due to the smallness of the anomalous dimen- 
sion i]*. (For the three-dimensional theory 77* ~ — 0.04, see section 9.) We note that for 
i]* = the fourth and higer derivatives of M*(p) vanish identically for p — ► 00. Numerical 
studies of eq. ( |4.10|) indicate that this is approximately true for small values of p as well 
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TPfl . Moreover, infinitesimally close to the phase transition the system spends arbitrarily 



long "time" \t\ = — ln(fc/A) near the fixed point, before deviating towards the phase with 
spontaneous symmetry breaking or the symmetric one. As a result the value of rj which 
is relevant for the dynamics is the fixed point value 77*. There should, therefore, exist an 
expansion in the value of 77 - the 77 expansion - where the derivatives u'"' for n > 4 give 
only small corrections ~ 77. For p — > 00 we observe the relations 

pu^\p)= -^-u"'{p) = -eu"'{p) (7.3) 

pV 5 )(p) = e{e+l)u"'{p) (7.4) 

and similar for higher derivatives. All derivatives with n > 4 are of order e ~ 77 
and combinations of the type p 2 u^ + pu^ are even suppressed by higher powers of e 
according to 

A truncation of ( f4.5|) should give reliable results at least for small values of 77*. Let us now 
try to exploit these properties also near the minimum of u and device a systematic trun- 
cation of eq. (|4.5|) . The evolution of u',u" is given by equations (|4.6[) , ( |4.7|) respectively. 
For the next two derivatives one finds 



dvl" 
~~dt 



(2d -6 + 3r])u"' + (d - 2 + r])pu {4) 
+ v d (N - \){2[u'fLl{u') - 3u"u"'L d (u') + u^L d {u')} 

+ v d {2(3u" + 2pu"'fL d (u' + 2pu") - 3{3u" + 2pu"')(hu"' + 2pu^)L d (u' + 2pu") 
+ (7u (4) + 2pu {5) )Li(u' + 2pu")) (7.6) 

(3d - 8 + 4t7)m (4) + (d - 2 + 7])pu {5) 

- v d (N - l){Q[u'fLi(u') - l2[u"] 2 u"'L d (u') + (Au"u w + 3K"] 2 ) L d 2 {u') - u^L d {u')} 

- v d {6(3u" + 2pu'"fL d {u' + 2pu") - \2{3u" + 2pu"') 2 (5u"' + 2pu^)L d (u' + 2pu") 
+ [A(3u" + 2pV")(7u (4) + 2pV 5) ) + 3{hu'" + 2pV 4) ) 2 ] L d (u' + 2pu") 

- (9« (5) + 2pu^)L d (u' + 2pu")). (7.7) 

As discussed earlier we can use these equations for a Taylor expansion of u around some 
fixed value. We choose p = po(k) if the minimum of u is away from the origin and p = 
otherwise. 

In the spontaneously broken regime (po{k) 7^ 0), the evolution of the minimum k of 
u(p) is given by eq. ( |4.8| ). We consider the couplings A = u"{k), u n = u^ 1 \k) which obey 
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the evolution equations 



du n du^ 

Hi 



+ u rM -. (7.8) 



du 3 



dt ^ dt 

The running of A is given by eq. ( |4.9|) and the next two couplings evolve according to 

& = (2d - 6 + 3 V )u 3 + v d (N - 1){2A 3 L^(0) - 3Xu 3 L d 2 (0)} 
+ v d {2{3\ + 2Ku 3 ) 3 L d 3 {2\n) - 3(3A + 2ku 3 )(5u 3 + 2kw 4 )L^(2Ak) 
+ [au, + 2ku 5 - <2 ^±\ L d (2\K)} (7.9) 

/3 4 = (3d - 8 + Ar])u A - v d (N - 1){6A 4 L^(0) - 12A 2 w 3 L^(0) + (4Aw 4 + 3u 2 3 )L 
- v d {6(3\ + 2Ku 3 ) A L d A (2\n) - 12(3A + 2ku 3 ) 2 (5u 3 + 2K Ui )L d (2\K) 



dui 
~dT 



+ 



4(3A + 2ku 3 )(7u 4 + 2ku 5 ) + 3(5u 3 + 2km 4 ) 2 L%(2\k 
2ku 3 u^\ ~ 



- [6u 5 + 2ku 6 - fH^j Li(2\K)}. (7.10) 



We observe that eqs. (|4.8|) , ( |4.9| ), (|7.9|) , ( |7.10|) in leading order in iV, i.e. keeping the 



contribution (n(ci — 2 + if) — d)u n and the term ~ (N — 1), can be solved level by level. 
In fact, /3 K , j3\ involve only k and A, (3 3 depends only on k, A, u 3 and /3 4 on k, A, u 3 , w 4 . 
In this limit no truncation is necessary for the determination of the infrared fixed point 
«*, A*, W3*, u 4 *. This behaviour generalizes to arbitrary n. The problem arises from the 
last bracket in the expressions for /3 3 , /3 4 which involves the couplings U5, ^6 which are 
undetermined at this stage. The same problem shifts to higher couplings if higher f3 n are 
computed. The function (3 n always involves the couplings u n+ i, u n+2 in addition to u m 
with m < n. At this place some truncation of the infinite system of coupled differential 
equations is required. We need some information on and uq. A first possibility is to 
put all u n with n > 4 to zero (truncation I). A second choice (truncation II) is suggested 
by eq. ( |7.5| ). We approximate 

u$ = — K^Ui, uq = 2k~ 2 m 4 for 0.01 < 2Ak < 100 

u 5 = u 6 = for 2Xk < 0.01 or 2Xk > 100. (7.11) 

As we shall see in section 9, at the fixed point 2k*A* ~ 1, independent of N. Through 
the truncation II we approximate w 5 , w 6 by expressions inspired by the analysis of the 
scaling solution, for the part of the evolution during which the system is near the fixed 
point. For the rest of the evolution we neglect them. Numerical values for the fixed point 
for the two different truncations are given in section 9, where all the numerical results are 
presented. 



22 



In the symmetric regime (po(k) = 0) we use the couplings m 2 = u'(0), A = u"(0), 
u n = u( n '(0). Their evolution equations read (in the symmetric regime z = 1, see eq. 
( |3.15D , and the L d , L d integrals coincide) 

dm 2 



dt 
d\ 
~dt 



Pm = (-2 + r/)m 2 + v d (N + 2)\L d {m 2 ) (7.12) 
P x = (d - 4 + 2 V )\ - v d (N + 8)\ 2 L d 2 (m 2 ) + v d (N + A)u^L d {m 2 ) (7.13) 



dU3 ~ fo = (2d-6 + 3r])u 3 + 2v d (N + 26)X 3 L d (m 2 ) - 3v d (N + U)\u 3 L%(m 2 ) 



dt 

du. 



+ v d (N + 6)u 4 L d (m 2 ) (7.14) 
p A = (3d - 8 + Ar])u 4 - 6v d (N + 80)A 4 L^(m 2 ) + 12v d (N + U)X 2 u 3 L d (m 2 



_4 _ q _ ioj o i a„\.. a.. I at i on\\4:Tdr^2\ , 10 „. I at i aa\\2„. 

dt 

- Av d (N + 2ti)\u 4 L d (m 2 ) - 3v d (N + 2A)ulL d (m 2 ) + v d (N + 8)M 5 L?(m 2 



(7.15) 

For the numerical work we shall neglect in the above equations, an approximation 
consistent with both truncations I and II. 

The last necessary ingredients for the numerical study of the evolution equations are 
the anomalous dimensions 77, rj, which we compute in the following section. 



8. Anomalous dimensions 

In this section we compute the anomalous dimensions 77 and fj (defined in eqs. ( |3.9|) ), 
as well as the evolution equation for the ratio z = Z^jZ^ (defined in eq. ( |3.15| )). The 
starting point is again the exact evolution equation ( |2.6|) , with Tk parametrized according 
to eq. ( 2.11 ). In order to obtain the evolution equation for (for N > 2), we consider a 



background field configuration with a small momentum dependence given in momentum 
space by 

0i(O) = 0, 5fa(Q) = 6<j>, (8.1) 

with 5(f) <C 0, and <p a = for a > 2. (We remind the reader that we are concentrating on 
real fields for which <fi(—Q) = 4>*(Q).) Inserting the above configuration into (|2.11|) we 
obtain 

ZM = Z t (p,Q> = 0) = ^jL 2l JEt^ {M = 0). (8.2) 



Similarly the evolution equation for Z^ is obtained from eq. ( |2.6| ) by considering a con- 
figuration 

0i(O) = 0, 5MQ) = ^ ( 8 -3) 

with 5(p <C (f>, and (f> a = for a > 2. Eq. ( |2.11| ) then gives 



Z M = Z t(p ,(? = 0) = l^^^-M = 0). (8.4) 
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For the evaluation of the r.h.s. of eq. 
action around the configurations (|8.1 



we need an expansion of the effective average 
^3|). This calculation has been done in ref. 114 . 



The results are summarized in appendix B, where the general expressions for rj and fj 
are also derived. In this section we limit our discussion to the approximation of uniform 
wave function renormalization introduced in section 4 (see eq. (|4.1| )). This approxima- 
tion corresponds to considering only the value of Z k (p), Z k (p) at the minimum po(k), 
while neglecting their derivatives. The relevant expressions are obtained by inserting the 
truncations 



Z'k(po) 



zl< 



into eqs. 



>fc(Po) = 

Z' k (p ) + Y k (p ) + Y k \po)po = 

Z' k \ Po ) + 2Y k \p ) + Y k "(p )p = (8.5) 

$B?2% ) and (|T09| ) (where the primes denote p-derivatives) . It is convenient to 
introduce the dimensionless quantity 

y = Z k 2 k d ~ 2 Y k (8.6) 

such that 

z = = 1 + Ky Q . (8.7) 

It should be noted that for N = 1 there is only one kinetic invariant. The wave function 
renormalization should be identified with Z k , while Y k can be set to zero. 

In the spontaneously broken regime one finds for the anomalous dimension r\ 



i 1 = --±K- 1 M*{2\K)-v d y 



L((2Xk) 



(8.8) 



and similarly for fj 

fj = v d (N - l)yo2 _1 
4v d 



L((0) + 2\nLi(0) 



d 



The functions M$ , 



— Z K 



(N - l)X 2 Z 2 k M^(0) + (3A + 2u z k) 2 Z{M%{2\k) 



(8.9) 



Mf n and M$ are defined in eqs. 



'4,0 



(|BJ6l) , ( gg|) and ( p^7|) re- 
8.9|) we are using rescaled arguments with the 



spectively. Notice that in eqs. 
conventions of eq. (|3.22|) . We shall need the function M| for the numerical calculations 
of the next section. Within the approximation of uniform wave function renormalization 
and with the notation introduced in section 3, it can be written as 




1 



^1 + r)y [z + r)y + w 




d 



dr 



dy 



( dr" 
- y [ 1+r + y dy- y 



y— r + rjr + r\y— 



dy 



8r" 

rjr + 2y- 



+ 



(1 + r)y (z + r)y + w 



[8.10) 
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(The dependence of M$ (and the other functions) on r] and z has been suppressed in 
the notation of eq. (|8.8| ) for simplicity.) This integral simplifies considerably if the 
contributions proportional to 77 are neglected and we set z — 1. With these simplifications 
it depends only on w and is discussed in the appendix C. The evolution equation for z 
can be expressed in terms of rj and fj 

dz d{Ky ) 



dt dt 



Hv-v)- (8-ii) 



For the scaling solution z is at its infrared stable fixed point z* (see discussion at the end of 
section 3). From eq. ( |8.11|) we conclude that the fixed point values r/* and fj* coincide. In 



other words, z adjusts itself so that the two anomalous dimensions become identical. This 
behaviour is expected to be a general feature, independent of the particular truncation 
used in this paper. 



In the symmetric regime, where po(k) = 0, we obtain from eq. (p.40|) 

V = V = -v d y Li{m 2 ). (8.12) 

At the transition between the spontaneously broken and symmetric regime (for k = 0, 
fh 2 = 0) this expression agrees with eq. ( |8.8| ) but not with eq. ( |8.9| ). The reason is that 
for small k the truncation Z' k (p ) = ceases to be consistent with a nonvanishing Y k (p ) 
(and Z' k (po) = 0). A reliable computation of fj for very small nonvanishing k should 
improve over the truncation of eq. ( |8.5| ). For similar reasons the present truncation does 
not guarantee z — > 1 for k — > in the spontaneously broken regime. Such a behaviour 
is required by continuity in a treatment without truncations, since z exactly equals one 
in the symmetric regime. We note that beyond the approximation Z' k (p ) = Z' k (p ) = 
the values of rj or fj do not necessarily coincide at the transition between the two 
regimes, since the definition of the anomalous dimensions involves additional terms from 
the k dependence of po(^) m the spontaneously broken regime. For very small k we may 
approximate fj in the spontaneously broken regime by neglecting the contribution from 
Z^ of eq. qgOg) and setting Z' k {p ) = 0, Y>{p )p = 0, U'"(p )p = 

fj = v d z- 1 y {2Li(2\K) + (N — 1)L?(0) + ^Lf + 2 (2Ak)}. (8.13) 

Finally we need the evolution equation for yo in the symmetric regime 

^= (d-2 + 2r ] )y + Z k 2 k d - 2 ^ (8.14) 



dY k 1 d 

hm — — 



dt p->o p dt 



Z k {p)-Z k {p) . (8.15) 



Here dZ k (p)/dt and dZ k (p)/dt can be extracted from eqs. (|B.7| ), ( |B.8| ) and ( |B.31| ), 
( |B.32| ) respectively, with W\ — Wi — U k (p) and Z k (p) = Z k (p). We observe that eq. 
( |8.15| ) is well defined since 

Hm ^Z fe (p) = l™J-Z k (p) = v.k^Z, 1 [NZ' k (0) + Y k (0)} L^rh 2 ). (8.16) 

P^O (Jt P^O (jt 
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Due to the presence of massive modes, quickly stops running in the symmetric regime 
and jjq approaches zero for d > 2. Since in the present paper the running in the symmetric 
regime plays only a secondary role, we shall use in this regime a somewhat stronger 
truncation by setting Yfc(O) = 0. This implies for the symmetric regime 

77 = 77 = 0. (8.17) 



9. Critical exponents of the three-dimensional theory 

In this section we use the formalism developed in the previous sections in order to study 
the phase transition of the three-dimensional A-component 4 theory. As we have already 
discussed in section 2, the effective average action interpolates between the classical action 
(for k = A) and the effective action (for k = 0). Therefore, our strategy is to define the 
classical theory at the ultraviolet cutoff scale and solve the evolution equation for the 
effective average action (in some appropriate truncation scheme), in order to determine 
the renormalized theory in the limit k — > 0. Subsequently, we can adjust the classical 
(bare) parameters of the theory in order to approach the fixed point in the critical region 
and obtain all the quantitative information for the phase transition. 

The evolution equations for the parameters of the effective average potential are given 
by expressions (|4.8| ), (|4.9|) , ( |7.9|) , (|7.10f) for the spontaneously broken regime, and ( |7.12|) 



- ( 17.151) for the symmetric regime. The derivatives of the potential higher than the fourth 
are either neglected (truncation I and symmetric regime) or estimated according to eq. 
( [7.1 1|) (truncation II). For the numerical calculation we consider only the anomalous 



dimension for the Goldstone modes rj. We have shown in section 8 that, for the scaling 
(fixed point) solution, the anomalous dimensions for the radial and Goldstone modes 
coincide: 77^ = 77^. Since in this work we are interested only in the critical region, it is a 
good approximation to drop the distinction between the two anomalous dimensions. For 
practical purposes we use eq. ( |8.8|) , keeping only the first term and neglecting the ones 
proportional to yo- It should be noted that this approximation is not valid for N = 1, 
where the only anomalous dimension is fj. It is expected, however, from continuity of our 
expressions in the parameter N, that eq. ( |8.8|) should give a good approximation to 77* 



even for N = 1. In the symmetric phase we set 77 = according to eq. ( |8.17[ ). For the 



numerical calculation we also drop the distinction between the integrals and as 



explained in section 4 after eq. ( |4.4| ). In this approximation these integrals are discussed 
in appendix A. Since they cannot be evaluated in closed form for a general argument, we 
perform numerical fits of L^w), which connect the expressions ( |A.4[) for w = with the 



asymptotic expressions ( |A.6| ). Another numerical fit is used for the integral M|(w) which 
is discussed in appendix C. 

A typical example of the numerical integration of the evolution equations is displayed 
in fig. 2, for N = 3. The truncation II is used, even though the parameters u^, U5, 
u e are not plotted. The evolution starts at k = A with the classical action, so that 
Z\ = 1, 113(A) = «4(A) = «s(A) = Uq(A) = 0. The quartic coupling is arbitrarily chosen 
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A(A) = A (A) /A = 0.1. Two values of k(A) are selected so that the system is very close 
to the two sides of the critical line. The two values are indistinguishable in the graph 
and are given by k(A) = K cr + 5k^, with K cr = [p (A)] cr /A ~ 0.107, \5ka\ <C K cr and Sk,\ 
positive or negative. With these initial conditions the system evolves towards the fixed 
point in the critical region and remains close to it for a long "time" \t\ = — ln(/c/A). The 
fixed point values of the different parameters are listed in the last row of table 1. By 
staying near the fixed point for several orders of magnitude in t, the system loses memory 
of the initial conditions at the cutoff. As a result the critical behaviour is indepenent of 
the detailed short distance physics and is determined by the fixed point, as long as the 
evolution starts sufficiently close to the critical line. In the final part of the evolution, as 
k — > 0, the theory settles down either in the phase with spontaneous symmetry breaking 
or the symmetric one, for 5k\ positive or negative respectively. 

a) For positive Sk\, the physical parameters of the theory approach their renormalized 
values in the limit k — > 0. Due to the presence of the mass for the radial mode, the 
anomalous dimension r\ becomes zero for k — > 0, and the wave function renormalization 
approaches a finite value Z = Z. We define 



POR = 


lim kn{k) = Z 1 p 

k^O 




\r = 


lim k\(k) = Z~ 2 \ Q 






limk 3 - n u n (k) = Z- n ul ) n \po). 


(9.1 



For N ^ 1 and k — > 0, the quantity X(k) approaches a finite value. As a result, due to 
the fluctuations of the Goldstone modes, the renormalized quartic coupling Xr is zero in 
the phase with spontaneous breaking and the theory is infrared free. This is not the case 
forN= 1, where A diverges in the infrared so that Xr remains non-zero, 
b) For negative 5ka, the quantity n{k) becomes zero at some non-zero k s . The evolution 
is continued with the equations for the symmetric regime, with m 2 {k s ) = and u n (k s ) 
assuming their values at the end of the running in the spontaneously broken regime. This 
part of the evolution is not displayed in fig. 2 for simplicity. Notice that the anomalous 
dimension approaches zero at k s and remains zero in the symmetric regime. For k — > 
the physical parameters of the theory approach their renormalized values 

m 2 R = lim k 2 m 2 {k) 

\r = lim k\{k) 

k—>0 

U nR = Yimk z - n u n {k). (9.2) 

Very close to the phase transition, rriR and \r are very small compared to A, with a fixed 
value for the ratio ttir/Xr. We shall return to this point later in this section. 

In table 1 the fixed point values of various parameters are listed, for N = 3 and various 
truncations of the evolution equations. For the first row only two evolution equations are 
considered for k and A, and the higher derivatives of the potential as well as the anomalous 

27 



dimension are set to zero. Then, the number of evolution equations is enlarged to three 
with the inclusion of u 3 . Subsequently 77 and w 4 are added (truncation I), and finally u 5 
and uq are approximated according to eq. ( |7.11| ) and are included in the system of five 



coupled evolution equations (truncation II). It becomes apparent from table 1, that the 
influence of the higher derivatives of the potential on the fixed point, and therefore on the 
critical dynamics, is small. In fact, for an efficient truncation it would suffice to consider 
the "relevant" first three p-derivatives of the potential and the anomalous dimension. The 
inclusion of higher derivatives gives small improvements in agreement with the analysis 
of section 7. In table 2 we list the fixed point values of the various parameters for several 
values of N. It is apparent that the anomalous dimension at the fixed point is always a 



small quantity, and that 2A*/t* ~ 1 for all N. This justifies eq. ( |7.11| ) for the truncation 
II. 

The critical exponents parametrize the singular behaviour of the theory as the phase 
transition is approached. A measure of the distance from the phase transition is the small 
difference 

Sk\ = k(A) — K cr (9-3) 

for any given value of A (A). If k(A) is interpreted as a function of temperature in a 
realistic three-dimensional model, the deviation 5k\ is proportional to the deviation from 
the critical temperature, = A R (T cr — T). The critical exponent (3 parametrizes the 
behaviour of the "unrenormalized" order parameter 

Po = Vanpo(k) = Z~ x p QR (9.4) 

according to 

Po ~ (^a) 2/? , (9.5) 

for 5k,\ — > + . Similarly, the critical exponent v describes the divergence near the phase 
transition of the correlation length, which is equal to the inverse of the renormalized 
mass rriR. As we have already mentioned, the theory is infrared free in the phase with 
spontaneous breaking for JV^ 1. This means that \ R vanishes and, since po R is fixed, the 
renormalized mass m R vanishes [] As a result, for Sk\ > the critical exponent v can be 
defined only for N = 1, according to 



m R = \j2\ R p 0R ~ (5k a ) u -. (9.6) 
For 5k\ < the theory is in the symmetric phase and, for any N, we can define 

m fl ~|«J« A r + . (9.7) 

(Here m R is given by m(k = 0).) The critical exponent 7 is related to the response 
of the expectation value <p a to an external magnetic field or source B a . The magnetic 
susceptibility 

d(b a , . 

X = (9-8) 



* We could avoid this problem by defining and m# by renormalized four and two point functions 
with non- vanishing external momenta. 
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evaluated for a vanishing magnetic field B a = 0, is non-analytic for 5ka — > 

X~N _1 - (9-9) 

The relation between cf) a and B a is encoded in the effective potential and reads for small 
values of (fi a and 5k a ^ 

B« = ®M = ^ = ^ 

X= rh' 2 . (9.10) 
Here the "unrenormalized" mass term rh 2 is given by 

m 2 = limm 2 (A;) = Zm 2 R . (9-H) 

k— >0 

For N — 1 and <5ka > we obtain 

m 2 = 2A p = PoR ~ (^a) 7 ~ . (9.12) 

In the symmetric phase (5ka < 0) we have for arbitrary iV 

rh 2 = Zm% ~ |5k a | 7+ . (9.13) 

For the numerical determination of the critical exponents we integrate the evolution 
equations with the truncation II, for several very small values of 5k\. The critical expo- 
nents are then obtained as the slope of the logarithm of the relevant quantity at k — 0, 
when plotted as a function of In \ 5k a \. More specifically, for the determination of (3 we 
compute the slope of In po as a function of In Sk\, and similarly for the other critical expo- 
nents. The behaviour of the exponents for decreasing \8k\\ is depicted in fig. 3 for N = 4. 
As Sk a — > they approach universal values. The deviation from the universal behaviour 
for larger values of 5k,\ depends on the details of the bare theory. Fig. 3 corresponds to 
a bare theory (defined at k — A) with N = 4, A (A) = 0.5, n cr ~ 0.120. We have verified 
that the asymptotic values of the critical exponents are universal by determining them 
for different A (A) and n cr . For the anomalous dimension the relevant value is the fixed 
point value r/*. The results of our calculation are summarized in table 3 for various N. 
For N = 1 the critical exponents v and 7 satisfy v + = and 7+ = 7_ with an accuracy 
of 0.1%. For N 7^ 1 we list the values for v + and 7+. The critical exponents satisfy the 
scaling laws 

P= (l + ?7*)^ (9-14) 
7= (2-7j> (9.15) 



with an accuracy of 0.1%. We mention also the scaling law of eq. ( |6.7|) , which gives the 
critical exponent 8 and was obtained in section 6. In table 3 we have included values for 
the critical exponents obtained through other methods (summed perturbation series in 
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three dimensions, e-expansion, lattice calculations and 1/iV-expansion), for comparison. 
We observe agreement at the 1-5 % level for the exponents (5, v, and 7. For the anomalous 
dimension 77* we observe satisfactory agreement of our results with the quoted values, even 
though 77* is a small quantity and most severely affected by the various approximations 
in our method. In fact, significant improvement for all the results is expected if the 
anomalous dimension for the radial mode fj is included. Preliminary results indicate 
that this requires a truncation that includes at least the first p-derivative of the wave 
function renormalizations. We postpone this discussion for a future publication |19| . 



An interesting point, which is connected to the resolution of the problem of infrared 
divergences of perturbation theory J13, IS[ |1^J, concerns the vanishing of the renormalized 



quartic coupling Xr at the phase transition. More specifically, the ratio Xr/tur takes an 
iV dependent value in the critical region. As a result, the behaviour of Xr in the critical 
region can be characerized by a critical exponent ( defined according to 

Xr ~ |^a| C (9.16) 

and 

C = v. (9-17) 

In table 4 we list our results for the ratio Xr/tur, as well as results obtained with other 
methods for comparison. For N ^ 1, due to the vanishing of Xr in the phase with 
spontaneous symmetry breaking, the ratio was calculated in the symmetric phase. For 
N — 1, the values obtained in the phase with spontaneous symmetry breaking and the 
symmetric one agree at the 0.1 % level. 



10. Conclusions 

The effective average action r*. results from the effective integration of fluctuations with 
characteristic momenta larger than a given infrared cutoff k. It is the appropriate quantity 
for the study of the physics at the scale k. It interpolates between the classical action 
S for k = A (A being the ultraviolet cutoff of the theory, much larger than any other 
physical scale) and the effective action T for k = 0. The dependence of on k is given 
by the exact evolution equation (|2.6| ) [l|]. No ultraviolet or infrared divergences appear in 
this evolution equation, and the formalism leads to an efficient treatment of the infrared 
problems which plague theories with massless modes in less than four dimensions. In 
this work we have studied the exact non-perturbative evolution equation for the effective 
potential. It is well suited for the study of critical phenomena since it can be cast into 
a scale independent form. In particular, a second order phase transition corresponds to 
fixed points for the various couplings appearing in the potential. They can be computed 
as the solution of a non-linear second order differential equation for the field dependence 
of the potential. We have given the explicit solution of this differential equation in the 
large N limit of the 0(iV)-symmetric scalar theory in three dimensions. It describes the 
fixed points for infinitely many dimensionless couplings. We also have solved the more 
general evolution equation away from the phase transition. 
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For a finite number of components a truncation to a finite number of couplings becomes 
necessary and the evolution equation has to be solved numerically. We have concentrated 
on the phase transition of the three-dimensional theory. Our results for the critical ex- 
ponents P, v, 7, 5 agree at the 1-5 % level with the results of the most sophisticated 
calculations through other methods (summed perturbation series in three dimensions, 
e-expansion, lattice calculations, 1/iV-expansion) . The anomalous dimension r\ is also 
well determined, even though significant improvement is expected with a more efficient 
treatment. We also emphasize that, apart from the precise determination of the criti- 
cal dynamics, the method of the effective average action permits the discussion of the 
non-universal behaviour of the theory, such as the behaviour away from the critical tem- 
perature for the finite temperature four-dimensional theory fL5 |. In this way permits the 
determination of the region of attraction of the possible infrared fixed points. The effective 
average action is not only useful for a study of fixed points. More generally it describes 
the relevant physics at a given lenght scale k~ l . In practice, this can be used through the 
identification of the infrared cutoff scale k with the appropriate physical infrared cutoff 
scale of the problem (such as the critical bubble scale for first order phase transitions or 
the Hubble parameter for inflationary cosmology). 

We conclude that the formalism of the effective average action provides intuitive un- 
derstanding, as well as computational power, for a large number of physical problems. 
Work in progress focuses on the equation of state for the 0(iV)-symmetric theory |19| , 
the further improvement of the precision for the calculation of the critical exponents, 
and the study of scalar theories with different symmetries and possible first order phase 
transitions J22|. The concept of the effective integration of degrees of freedom through 



averaging has been developed for fermionic and gauge fields |23], and an exact evolution 



equation has been formulated for gauge theories ||24j| . Work in this direction will lead 
to the study of more complicated systems, with exciting possibilities for new physical 
behaviour. 
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Appendix A: The integrals L 



In this appendix we disuss the integrals L d (w) = L d (w,0), where L d (w,rj) are defined in 
eq. O with L d (w,r},z) given by eqs. With the choice fl27|) - (Q for i4(g), 

they read 

-i 



+1 



exp(-y) 



[1 -exp(-y)] 2 
exp(-y) 



2n/ 2 
[1 -exp(-y)J 



1 - exp(-y) 

y 

1 - exp(-y) 



+ w 



1 -(n+l) 



+ W 



for 7i > l(A.l) 



They obey the relations 



nL d n+1 {w) 



for n > 1. 



(A.2) 



For u> > —1, the integrals |L^(u>)| are finite, monotonically decreasing functions of w |. 
We define 

L d (0) = -2l d n (A.3) 
and give the values of l d for the first five values of n 



'o 



r 

n 



jd 
'3 



I' 1 



CC= + i V ( d - + i 



'd 



2 (i-2-i)rg-i) 
3 (i_ 2 3-f + 3 2 -f) r g-2^ 

4 (l - 3 x 2 3 -i + 3 4 ^ - 4M) r g " 3 ) • 



We define the threshold functions s d (w) through 



Lt(w) 



<s d (w). 



(A.4) 



(A.5) 



They are decreasing functions of w, and 1 > s d (w) > for < w < oo. The first two 
terms of the asymptotic expansion of L^iw) for large arguments w — > oo read 



L d (w) 



' d 



' d 



2C hr + 1 T - + 2 \w 



r.-l 



C o + i +CU + 2 



§ The pole structure of L d n {w) is not relevant for this work. For a discussion see ref. [pi 
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L d Jw) 



+ n(n + l) 



C {l + X ) + C it + 2 )] 1 it + 3 ) *" (B+2) - fOT 72 " L (A - 6) 



Appendix B: Wave function renormalization 

In this appendix we derive general expressions for the anomalous dimensions rj, fj. We 
consider an effective action of the form of eq. (|2.11j ) . This is no longer the most general 
ansatz contributing to the wave function renormalization and our equations are therefore 
not exact, in contrast to the ones for the effective average potential. In addition, we 
neglect the q 2 dependence of the couplings Z k (p,q 2 ). This truncation corresponds to the 
most general action containing no more than two derivatives. For notational simplicity, 
we denote in this appendix U k (p), Z k (p), Z k (p), Y k (p) by U, Z, Z, Y respectively. Primes 
on U, Z, Z, Y denote derivatives with respect to p. The wave function renormalizations 
Z k , Z k , Y k are defined at the minimum of the potential po(k) according to eq. ( |3.8|) and 
the anomalous dimensions according to eqs. ( |3.9| ). Also primes on Z k , Z k , Y k denote 
p-derivatives of Z k (p), Z k (p), Y k (p) at p = p . We point out that for N = 1 one should 
set Y = and Z is the relevant wave function renormalization. 

For a calculation of the evolution of Z (which is connected to the anomalous dimension 
rj) we need an expansion of the effective average action T k around the configuration given 
by eq. ( |8.1|) . This calculation can be found in section 5 of ref. |^4| . Here we present the 

(2) 

result for the second functional derivative T k . It can be written as 



(r k 2) )ab{q,q') = (Fo)ab(q,q') + (^i)ab(q,q') + {T 2 ) a b{q,q'), 

where we have separated the terms containing zero, one and two powers of 
term reads 

(r U(<7, q') = [(Z5 ab + P Y5 al 5 bl ) q 2 + M 2 ab ] 8{q - q'), 

where 

M 2 ab = U'5 ab + 2pU"5 al 5 bl . 
The second has only an off diagonal contribution 



(B.l) 

The first 

(B.2) 
(B.3) 



with 



(T 1 ) 22 (q,q')=0, (T 1 ) ab (q,q') = for a, b > 2 
f i(g, Q)5(q ~q'-Q) + Tl(q, -Q)S(q - q' + Q) 

(rDi2(g',g) 
l. 



(ri)u(q,q') 

(ri)i 2 (g,g'; 
(ri) 2 i(g,g'; 

f x (g, Q) = <t>5cf>{U" + "-Yq 2 - Z\q - Q)Q) 



The diagonal elements of the third term read 



:b.4) 



(r 2 ) aa (g, q) = 5<P*5<P{U" + 2U"5 a2 + 2U'"p5 al 

+ q 2 (Z' + Y5 a2 + Y'p8 al ) + Q 2 (Z' + Y5 a2 + 2Z"p5 al )}. (B.5) 
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Subsequently one has to evaluate the trace on the r.h.s. of eq. Q2.6| ) and extract the piece 
proportional to Q 2 {5<p*84>), which gives the evolution of Z according to eq. (|8.2| ). The 



details for the calculation of the trace are given in ref. ||14|| . For the scale dependence of 
Z (= Zk(p)) at fixed p one obtains 



dZ_ 
~dt 



Of 



+ 



dZ® 



Of 



:b.6) 



with 



Of 
dZ<® 



dt 



v d [(N - l)Z' + Y]Z k k L d {w\) + v d (Z' + 2Z" p) Z k l k d ~ 2 L d (w 2 ) (B.7) 

4v d [U"} 2 pk d ~ e Q d 2 j(w u w 2 ) + Av d YU"pk d -*Q d 2 j(w l7 w 2 ) 
+ v d Y 2 pk d ~ 2 Q d 2 j( Wl ,w 2 ) - 8v d Z'U"pk d ~ 4 Ll 1 (w 1 ,w 2 ) 



4t> 
~d 



d [Zfpk d ~ 2 L d+2 (w u w 2 ) - Av d Z'Ypk d - 2 L d+2 (w u w 2 



+ ^Z'W'pk^N^iw^w,) + ^Z'Ypk d ~ 2 N^ 2 ( Wl ,w 2 ). 
a a 

Here we have used the variables 

wi = U' 

w 2 = U' + 2pU". 

It is also convenient to define the shorthands 

P = Zx + R k (x) 
P = Zx + R k (x) 

with Rk(q) given by eq. ( |2.2|) and x = q 2 . The integrals L d , L d are given by 

d, d 

in ^ 



(B.i 



CB.9) 



(B.10) 



POO , fj 

L d n (w)= k 2n ~ d Zl \ dxx*- 1 — (P + w)~ n 
Jo at 



L d Jw) 



POO , 3 ~ 

k 2n-d z n I dxX^ — iP + w)- 71 , 







dt 



[B.ll) 



where we emphasize that d/dt acts on the r.h.s. only on R k (x). In the approximation 
of uniform wave function renormalization (when the p dependence of Z, Z is neglected) 
these integrals are discussed in the beginning of section 4 and in appendix A. Notice the 
connection between the arguments w and w: w = w/Z k k 2 . The dimensionless integrals 

L d ni , n2 (w u w 2) = k 2{ni+n2) - d J o rfxx2- 1 -{(P + Wl )"" 1 (^ + ^2)-" 2 } (B.12) 
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are generalizations of the integrals L^, L d with 

L%( Wl ) = Z k n L d n ( Wl ) 

L d 0jn (w 2 ) = Z^L d n {w 2 ). (B.13) 
Similarly, the integrals N d n2 (wi,w 2 ) are defined by 

roo d 3 r f) P ~ -i 

N* una ( Wl ,W 2 ) = k *(m+n,-l)-d dxx - 2 _{_ { p + Wl) -n 1{p + W2) -n 2 y (B>M) 

Finally, the integrals 

Qn"n 2 (Wl, W 2 ) = ni+n2 ~ a) ~ d / dxx?- 1+a 



d_ 

dt' 



OP 2xd 2 P Ax fdP\ 
+ 1 1 



Ax (dpy , n > i 



dx d dx 2 

are related by partial integration to other integrals 



(P + Wi)- ni (P + w 2 )- na f&.15) 



, a 3 f fdP' 
dXX ~ 2 dt\ldx-. 



(P + w 1 )- ni (P + w 2 ) 



:b.i6) 



through 



+ ^ P YN^ +1 ( Wl ,w 2 )-^N% 



2a-2, 
n 2 ' 



Wi, W 2 . 



(B.17) 



The evolution equation for Z depends, in this approximation, on [/', U", Z, Z', Z" and 
Y . There is also an 77 dependence resulting from the t-derivative acting on Z)~ in Rk{q)- 
We define the anomalous dimension 77 by the variation of In Z^ = In Z/,(p (k)) (see eqs. 
[8|), ( |3.9| )). In the spontaneously broken regime, where po(k) ^ 0, this implies 



0t 



+ 



Po 



dt 



- Z^Z' k 8. 



fB.li 



PO/ 



Here 5 = dpo(k)/dt describes the k dependence of the minimum of the potential 



5 = - v d Z^k d ~ 2 {(N - l)Lf(O) + (3 + L((2X Po )} 
- v d Z k l k d {{N - 2 (0) + &f 2 (2A Po )}, 



(B.19) 



with 



A = ^( Po ), U n = ul n) (j>o) for n>3 
wi(po) = 0, u>2(po) = 2Ap - 



^B.20) 
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Combining eqs. QB77D , ( gg ), ( [BTjg ), (gTTgJ) we obtain 

77 = - ^- 2 Z fe - 2 {[(iV - l)^ + Y k ] Lf(0) + + 2p Q Z ! l)L d 1 (2\p Q )} 



,2\2 



^ 1 A%Wa^(0,2^) + ^Mg 2 (0,2A Po ) +^ (^-1 Mg 4 (0, 2A Po ) 



^^^^{^(O^Apo) + ^A$ 2 (0,2Ap ) + £ (^V^t 4 (0,2V ) 



- ^Z k \2Z' k - Y.yXpok'-^N^Xpo) + l^iV 2 d j 2 (0,2Apo)} 

+ SvdZ^Z'^pok^Lf^O, 2Ap ) + ^Z^Z' k (Z' k + dY k )p k d - 2 L d + 2 (0, 2A Po ) 

- Z^Ztf. (B.21) 
It is convenient to bring rj into the form 

rj= - ^Z k l p^k d ~ 2 {M d fi (0) -2M? ;1 (0,2Apo) +M d 2 (2Ap )} 

- ^Z fc - 2 n£; d - 2 {L?(2Vo) - ^ fc - 1 n.po^ +2 (2Ap )} 

- v d Z^Z' k k d ~ 2 ^{N - l)Z k x L d M - ^(0, 2A Po ) 

+ ^ (5 + Lf(2Ap ) - ^^Po^ 2 (0, 2Ap )| 

- Z^Z' k 8. (B.22) 
For this purpose we have employed the following identities 



k 2 



w 2 - Wi 



(dKtft) - <-i,n 2 K^2) - YpN d + 2 2 ( Wl ,w 2 )} (B.23) 



<o(0) = ^^ Z ^ 1)L n-M for n > 1 (B.24) 

<,nH = ^^^ (n_1) ^iH-^^n +2 H for n>l (B.25) 
M^ n2 (w 1 ,w 2 ) = 

: ^-^{M d i , n2 _ 1 (w u w 2 ) - M d ni ^ n2 { Wl) w 2 ) - YpM d + 2 m (w u w 2 )}. (B.26) 
We also define the function 



M 2 » = M d fi (0) - 2M d 1 (0,w) + M d 2 (w) 
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d d f dP 

•2 i 



k 2 d dxx2 n . 

Jo dt[dx \P P + w 



1 



1 



^B.27) 



For the calculation of the evolution of Z (which is connected to fj) we need an expansion 
of Tfc around the configuration given by eq. ( |8.3|) . The second functional derivative 1^ 



is again given by eq. (B . 1| ) with T given by eqs. ( B.2|) , ( B.3|) . The other two terms now 
read 

(FiUte, q') = (f i) a6 (g, Q)6(q -q'-Q) + (H Wg, -Q)S(q - q' + Q) 
(T 1 ) ba (q,q') = (nU(q',q) 

1 



with (Ti) ab (q,Q) 



U"+(q 2 -qQ)Z' + -Q 2 Y 



6 n 



+ 



1. 



2U" + 2U"'p + (g 2 - qQ)(Y + Y'p) + Q 2 (Z' + -Y + Y'p) 



<5 al <5 61 (B.28) 



and 



(r 2 ) a6 (g, q) = 6<p*6<p{(U" + 2U"'p)S ab + (2U" + 10U'"p + AU^p 2 )S al 5 bl 
+ (Z' + Y'p)Q 2 5 ab + (Y + 2Z"p + AY'p + 2Y" p 2 )Q 2 6 al 6 bl 
+ (Z' + 2Z"p)q 2 5 ab + (Y + 5Y'p + 2YV )g 2 5 al 5 w }. (B.29) 



For the calculation of the trace we refer the reader again to ref. |TJ| and we give the result 
for 



dZ_ 
dt 



dt 



+ 



dt 



[B.30) 



One finds (with wi, w 2 given by eqs. (B.9Q ) 



ot 



dzw 



v d (N-l)(Z' + Y'p)Z^k d - 2 Lf(w 1 ) 



+ v d {Z' + 2Z"p + Y + 5Y'p + 2Y" p 2 )Z^k d - 2 L{{w 2 ) 



(B.31) 



dt 



2v d (N - l)[U"] 2 pk d ~ 6 Q d 3 '° ( Wl ) + 4v d (N - W'Z'pk'-'QMw!) 



+ 2v d (N - l)[Zfpk d ~ 2 Qf 2 ( Wl ) + 2v d (3U" + 2U>" p) 2 pk d ^Q d ^{w 2 ) 

+ 4v d (Z' + Y + Y'p)(3U" + 2U"'p)pk d - A Q%(w 2 ) + 2v d (Z' + Y + Y' p) 2 pk d ~ 2 Q%{w 2 ) 

+ *£(N - l)Z'U"pk d - 4 N d ( Wl ) + ^(N - l)[Zfpk d - 2 N d + 2 ( Wl ) 



+ ^(Z' + Y + Y'p)(3U" + 2U">p)pk d -*N d fi (w 2 ) + ^{Z' + Y + Y' p) 2 pk d ~ 2 N d + 2 (w 2 ) 
- 2v d (N - l)YU"Z?pk d -*L*(wi) - 2v d (N - 1) [z'Y + ^[Z'] 2 ^ Z k 2 pk d - 2 L d 2 +2 ( Wl ) 
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Av d (Z' + Y + Y'p)(3U" + 2U"'p)Z^pk d - 4 Lj(w 2 



-2(2 + -) v d (Z' + Y + Y'p) 2 Z k 2 pk d - 2 L a 2 ■ - ( <r 2 ) 



Here we use the definitions 



(B.32) 



K un2 (w 1 ,w 2 ) = N* un2 (w u w 2 ;P^ P) = N^ ni (w 2 ,w 1 ) + YpL d + 2 ni (w 2 ,w 1 ) (B.33) 



and 

%°n>i. w 2 ) = Q^>i, w 2 ; P <-> P). (B.34) 

The integrals Q^ n „ are related to M% n2 by an equation analogous to eq. (B.17| ). The 
relation between M d „„ and M d „„ reads 

M^ n2 (w u w 2 ) = M d iin2 ( Wl ,w 2 ;P^ P) = 

M d n2m {w 2)Wl ) + 2YpN d 2tni (w 2 ,w 1 ) + Y 2 p 2 L d + 2 ni (w 2 ,w 1 ).(B.35) 

We also observe a general relation for the integrals N% n2 {wi, w 2 ) 

(m - l)N d i n2 { Wl ,w 2 ) = ~L d i _ ln2 (w 1 ,w 2 ) -n 2 iV"4 +lini _ 1 (w 2 ,w 1 ) = 

T^t-i^K^) - n 2 ypL^ lin2+1 (wi,«; 2 ) - n 2 N d 1 _ lin2+1 (w 1 , w 2 ).(B.36) 



From this, or directly from partial integration of eq. (|B.14j ), follow the identities 





» = 


d L d 
4^2,0 




K 


» = 


4^0,2 






» = 


4^0,2 


'w)-YpL d f(w) 


^0,3 


» = 


d L d 


[w) + YpL d f{w) 



^B.37) 



Using the various relations among the integrals we can express dZ/dt in terms of M d , 
M d , L d and L d only We conclude that the evolution equation for Z (= Z k (p)) depends 
on U', U", U'", Z, Z', Z", Y, Y' and Y" in our approximation. There is also an rj 
dependence resulting from the t-derivative acting on Z^ in Rk(q)- 

The anomalous dimension for the radial mode is defined in the spontaneously broken 
regime by 



Of 



+ 



Of 



B.38 



pa , 



38 



We find 

fj = - v d Z^Z^k d - 2 {(N - l)(Z' k + n'Po)^(O) + (Z' k + 2Z'lp Q + Y k + 5Y^p + 2Y k "pl)L d 1 (2\p )} 
- ^Z k l Pok d -'{(N - l)A 2 M 4 d (0) + 2(N - l)Z>\k 2 M d + 2 (0) + (N - l)[Z' k ] 2 k±M d +\0) 
+ (3A + 2U 3 p ) 2 Ml (2Xp ) + 2(4 + Y k + F fc 'p )(3A + 2U 3 p )k 2 M^ 2 (2Xp ) 
+ {Z> k + Y k + Yip ) 2 k i Mi+\2\p )} 

+ 2v d Z k 2 Z k l p k d ~ 4 S [ (N - l)Y k XL d 2 (0) + (N - 1) (zjft + ~[4] 2 ) fc 2 Lf 2 (0) 

+ 2(4 + n + Y k ' Po )(3X + 2U 3 p )L d (2Xp ) + (2 + i) (4 + K fc + F fc ' P o) 2 A; 2 Lf 2 (2Apo) | 
-4- 1 (4 + n + F fe V )5. (B.39) 

Finally, we need the anomalous dimensions in the symmetric regime where po(k) = 0. 
In this regime Z k and coincide. The contribution from dZ^ /dt or dZ^ /dt vanishes 
identically for p = 0, since there is no cubic vertex. One finds from eq. ([B.7| ) or eq. 

(HI 



77 = 77 



dlnZ fc (0) dlnZ fc (0) 



dt 



dt 



-v d Z k 2 (NZ' k (0) + Y k (0))k d - 2 L d (m 2 ) (B.40) 



with m 2 = U' k (0). 



Appendix C: The integral 

In this appendix we disuss the integral M d {w) = M d (w, 0, 1), where M d (w, rj, z) is defined 
in eq. (|8.10| ). Since z — 1, it can be written as 

M d (w) = w 2 Z 2 M d 2 (0,w) } (C.l) 

where M d x n2 (u>i, w 2 ) is defined in eq. ( B.16|) . With the above approximations (uniform 



wave function renormalization, rj = 0, z = 1) and the choice ( |2.2|) - ( |2.3|) for R k (q), this 
integral reads 



M 2 d 2 (0,w) 



g-2 



~i2 



(1 + r) 2 [(1 + r)y + w] 



2y— + 2 (y—) r - 2y 2 [l + r + y— ) — 



+ 



(1 + r)y (1 + r)y + w 

(C.2) 
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with 



r(y) = eM ~ y) 
{V> 1 - exp(-y) 



The function M| 2 (0,w) is monotonically decreasing for positive w. We define 

Mj 2 (0,0) = -2Z t :M 

and find 



mt = 2T --1 



d ( d 
1-2-3 1 + - 



We can define the function t 2)2 (u») through 

M*(0,w) = -2Z^ 2 mttUw). 



(C.3) 

(C.4) 
(C.5) 

(C.6) 



It is a decreasing function of w, and 1 > t 22 (u>) > for < w < oo. The first term of 
the asymptotic expansion of M| 2 (0,w) for w — > oo reads 



M 2 d 2 (0,*)=4Z fc - 2 r 



1 -d + 



d(d + 2) 



C 



'd\ (2-d)d (d 



8 



For d = 2 the leading term takes the particularly simple form 

Ml 2 (0,w) = -2Z^ 2 w- 2 ... 

and 

lim M 2 (w) = -2. 



C(o + lj|*" 2 - (C7) 

(C.8) 
(C.9) 
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Tables 



ft* 


A* 


u 3 * 


V* 


«4* 


u 5 * 


«6* 


6.57 x 1(T 2 


11.5 












8.01 x 10~ 2 


7.27 


52.8 










7.86 x 10~ 2 


6.64 


42.0 


3.58 x 10~ 2 








7.75 x 10~ 2 


6.94 


43.5 


3.77 x 10~ 2 


95.7 






7.71 x 10~ 2 


7.03 


43.4 


3.83 x 10~ 2 


111 


-1.43 x 10 3 


3.72 x 10 4 



Table 1: Fixed point values for the minimum of u, its derivatives and the anomalous 
dimension, for various truncations of the infinite system of evolution equations. N = 3. 



N 




A* 


u 3 * 


«4* 


V* 


2ft* A* 


1 


4.05 x 10~ 2 


9.25 


87.4 


598 


4.49 x 10~ 2 


0.749 


2 


5.81 x 10~ 2 


8.08 


61.5 


249 


4.23 x 10~ 2 


0.939 


3 


7.71 x 10~ 2 


7.03 


43.4 


111 


3.83 x 10~ 2 


1.08 


4 


9.71 x 10~ 2 


6.15 


31.3 


56.6 


3.42 x 10~ 2 


1.20 


10 


2.26 x lO" 1 


3.26 


7.60 


6.01 


1.87 x 10- 2 


1.48 


20 


4.49 x lO" 1 


1.77 


2.14 


0.944 


1.02 x 10~ 2 


1.59 


100 


2.24 


0.375 


9.32 x 10~ 2 


9.05 x ID" 3 


2.17 x 10-3 


1.68 



Table 2: Fixed point values for the minimum of u, its derivatives and the anomalous 
dimension, for various values of N, with the truncation II (eq. (|7.11j )). 
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N 


P 


V 


7 


V* 


1 


0.333 


0.3250(15) a 
0.3270(15)° 
0.312(5) c 


0.638 


0.6300(15) a 
0.6310(15)° 
0.6305(15) c 


1.247 


1.2405(15) a 
1.2390(25)° 
1.2385(25) c 


0.045 


0.032(3) a 
0.0375(25)° 


2 


0.365 


0.3455(20) a 
0.3485(35)° 


0.700 


0.6695(20) a 
0.671(5)° 
0.672(7) c 


1.371 


1.316(25) a 
1.315(7) 6 
1.33(2) c 


0.042 


0.033(4) a 
0.040(3)° 


3 


0.390 


0.3645(25) a 
0.368(4)° 
0.38(3) c 


0.752 


0.705(3) a 
0.710(7)° 
0.715(20) c 


1.474 


1.386(4) a 
1.39(1)° 
1.38(2) c 


0.038 


0.033(4) a 
0.040(3)° 


4 


0.409 




0.791 




1.556 




0.034 




10 


0.461 


0.449 d 


0.906 


0.877 d 


1.795 


1.732 d 


0.019 


0.025 d 


20 


0.481 


0.477 d 


0.952 


0.942 d 


1.895 


1.872 d 


0.010 


0.013 d 


100 


0.497 


0A96 d 


0.992 


0.989 d 


1.981 


1.975 d 


0.002 


0.003 d 



Table 3: Critical exponents of the three-dimensional theory for various values of N. For 



comparison we list results obtained with other methods as summarized in |L2| and |21 

a) From summed perturbation series in fixed dimension 3 at six-loop order. 

b) From the e-expansion at order e 5 . 

c) From lattice calculations. 

d) From the 1 /iV-expansion at order 1/N 2 . 
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N 


1 


2 


3 


4 


10 


20 


100 


Aa. 

m R 


9.6 


7.9" 


8.4 


7.1 a 


7.4 


6.4 a 


6.6 




3.8 


5.0 6 


2.1 


2.5 6 


0.49 


0.50 b 



Table 4: The ratio Xr/ttir for the three-dimensional theory in the critical region, for 
various values of N. For comparison we list results obtained with other methods: 

a) From summed perturbation series in fixed dimension 3 [0, [12| . 

b) From large N calculations |T7L [12 . 



Figures 

Fig. 1 Graphical representation of the exact evolution equation for the effective average 
action T^. 

Fig. 2 Solution of the evolution equations for n(k), X(k), u^(k) and rj(k) in the critical 
region, for initial values slightly above and below the critical line. The fixed point 
for the second order phase transition is displayed, as well as the final running towards 
the phase with spontaneous symmetry breaking or the symmetric one. N = 3. 

Fig. 3 The critical exponents as the phase transition is appoached. N = A. 
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